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Abstract 

X-ray polarimetry has the potential to make key-contributions to our un- 
derstanding of galactic compact objects like binary black hole systems and 
neutron stars, and extragalactic objects like active galactic nuclei, blazars, 
and Gamma Ray Bursts. Furthermore, several particle astrophysics topics 
can be addressed including uniquely sensitive tests of Lorentz invariance. 
In the energy range from 10 keV to several MeV, Compton polarimeters 
achieve the best performance. In this paper we evaluate the benefit that 
comes from using the azimuthal and polar angles of the Compton scattered 
photons in the analysis, rather than using the azimuthal scattering angles 
alone. We study the case of an ideal Compton polarimeter and show that 
a Maximum Likelihood analysis which uses the two scattering angles low- 
ers the Minimum Detectable Polarization (MDP) by ~20% compared to a 
standard analysis based on the azimuthal scattering angles alone. The ac- 
curacies with which the polarization fraction and the polarization direction 
can be measured improve by a similar amount. We conclude by discussing 
potential applications of Maximum Likelihood analysis methods for various 
polarimeter experiments. 
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1. Introduction 



Some of the most interesting astrophysical objects are too small to be 
spatially resolved by present-day telescopes: the Schwarzschild radius of a 
galactic stellar mass black hole of Mi = 10 Mq solar masses at a distance 
of c?3 kpc corresponds to an ang ular extent of 2x10-"^ Mi/rfg /iarcsec; for 
an extragalactic supermassive black hole of Mg = 10^ M© solar masses at a 
distance of 20 rfg Mpc, the Schwarzschild radius corresponds to an angular ex- 
tent of Mg / (ig /iarcsec; a time of At = 1000 sec (in the observer frame) after 
the explosion of a Gamma Ray Bursts (GRB) at redshift z = 1, a GRB jet 
has an angular extent of 1.2 x 10"^° iV/lQ^f (At/1000 sec) (sin (^)/sin (1°)) 
/iarcsec if we assume (somewhat unrealistically) a constant jet expansion ve- 
locity with a bulk Lorentz factor F along an angle 9 to the line of sight. The 
numbers demonstrate the difficulties to image these objects. In the absence 
of imaging information, X-ray spectropolarimetric observations can deliver 
the additional information required to decide between models which cannot 
be distinguished based on spectroscopic data alone 

An X-ray polarimeter can deliver for each time interval not only an energy 
spectrum of the X-ray flux F{E), but also the energy spectra of the X-ray 
polarization degree ao(-E) and the X-ray polarization direction (j)Q[E). X-ray 
polarimetry will make a big leap when NASA will launch the Gravity and 



Extreme Magnetism SMEX [GEMS] mission in 2014 |lO|. G^M^ will use 



two Wo Iter type mirror assemblies to focus 2-10 keV X-rays onto two photo- 
effect polarimeters. Each polarimeter is made of a photoelectron-tracking 
time projection chamber. A student experiment will give additional polar- 
ization sensitivity at 0.5 keV. 

In this paper, we discuss the analysis of data from polarimeters in the 
> 10 keV to several MeV energy range based on the detection of Compton- 
scattered photons. Compton polarimeters use the fact that photons scatter 
preferentially perpendicular to the orientation of the electric field vector of 
a polarized X-ray beam. Several past and present satellites were equipped 
with detector configurations which - although not designed for Compton po- 
larimetry - were used to constrain the polarization of the X-ray emission 
from cosmic sources based on the analysis of Compton events (see 0, a and 
references therein). Recently, two experiments on board of the INTEGRAL 
(International Gamma-Ray Astrophysics Laboratory) satellite were used to 
measure the polarization of the 100 keV to 1 MeV emission from the Crab 
Pulsar and the Crab Nebula. Events scattered between different Ge detectors 
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of the SPI (SPectrometer on INTEGRAL) detector were used to measure a 
polarization degree of 46%±10% and a polarization direction aligned with 
the orientation of the X-ray jet [lH. The analysis of events scattered be- 
tween different CdTe detectors of the IBIS (Imager on-Board INTEGRAL 
Satellite) detector confirmed the large polarization fraction of the hard X- 
rays ll2\. In a similar fashion, Compton polarimetry was used with the Ge 
detectors onboard of the RHESSI experiment to constrain the polarization 
of solar flares in the >200 keV energy range 13|, [ij]. A report of polar- 



ized 7-ray emission from a Gamma-Ray Bursts remained controversial after 



16, 17 . The 



several independent analyses did not confirm the result [15|, 
soft gamma-ray telescope on board of the Japanese/US ASTRO-H mission 
(launch foreseen in 2013) will be able to measure the polarization of hard 
X-rays by detecting Compton scattered photons with a combination of Si 
pad detectors and Cadmium Telluride (CdTe) pixel detectors [18l] . 

A considerable number of dedicated Compton X -ray p olarimeters have 
been described in the literature, see for example H B HQ H H H i, [isj 
and references therein. Whereas some polarimeters are designed to measure 
exclusively the azimuthal scattering angles (e.g. 



20I . |23[ ) , other experiments 



also give information about the polar scattering angle on event to event basis 



(e.g. [18|, I21I, |22| ) . The analyses used so far for both types of experiments are 
based on fitting a constant offset plus a sine curve to the azimuthal scattering 
angle distributions. In this paper, we explore for the first time the sensitivity 
gain that can be achieved when using both, the azimuthal and the polar 
scattering angles in the analysis with the help of a Maximum Likelihood 
analysis method. In Sects. |2] and [3] we motivate the use of the polar and 
the azimuthal scattering angles based on the properties of the Klein-Nishina 
cross section and describe the standard analysis method and the Maximum 
Likelihood analysis method. In Sect. H] we discuss simulations of an ideal 
Compton polarimeter and evaluate the sensitivity gain achieved with the 
Maximum Likelihood analysis method. In Sect. |5]we summarize the results. 



2. Compton Scatterings: Basics 

Compton scatterings are governed by Compton's formula and by the 
Klein-Nishina cross section. Compton's formula reads: 



AA 
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where h is Planck's constant, mg is the electron mass, c is the speed of light, 
and 6 the polar scattering angle. The Klein-Nishina cross section is given by 
(see e.g. [26|): 
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with ro the classical electron radius, kg and ki the wave-vectors before and 
after scattering, 9 the scattering angle (the angle between ko and ki), and rj 
the azimuthal scattering angle, i.e. the angle between the electric vector of 
the incident photon and the scattering plane. Averaging over the azimuthal 
scattering angles gives the cross section for a scattering of a polarized photon 
with a polar scattering angle 6 irrespective of r/: 
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This cross section equals the Klein-Nishina cross section for unpolarized in- 
cident photons. 

As linearly polarized X-rays are preferentially scattered perpendicular to 
the orientation of the electric field vector, the probability density function of 
the azimuthal scattering angle ri exhibits a modulation of the form: 



P{r]) = - (1 + 6cos 
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The modulation amplitude b of the azimuthal scattering angle distribution 
can be computed from the Klein-Nishina cross section: 
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(S.=./2+(S).=o ko/k, + k,/ko 

Fig. [1] shows the modulation amplitude as function of 6 for several initial 
photon energies. The panel shows that events with scattering angles close to 
tt/2 exhibit a larger modulation amplitude. As a consequence, those events 
have a higher information content than events with scattering angles close 
to and close to vr. 



3. Analysis Methods 

In the standard analysis the azimuthal scattering angles 
togrammed. The polarization is determined by fitting the model 

ra(ao,0o;0) = n(l -haoyw cos [2(0-00)]) 



are his- 
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Figure 1: The modulation amplitude b of the azimuthal scattering angle from Compton 
scatterings as function of the polar scattering angle for photon energies of 10 keV, 100 
keV, 500 keV, and 1 MeV (from top to bottom). 



with a x^"fit to the histogram. The value n is the mean number of entries 
in each bin of the histogram: n = Uq/N, if tiq events were recorded and the 
histogram has N bins. The fitting parameter oq is the fractional polarization 
of the beam and equals 1 for a 100% polarized beam and for an unpolarized 
beam. The parameter n is the "modulation factor" which gives the modu- 
lation amplitude for a 100% polarized beam. In general, n depends on the 
physics of Compton scatterings and on the design of the polarimeter. In the 
following we will consider an "ideal polarimeter" which detects all Compton 
scattered photons and gives the scattering angles with very high precision. In 
this particular case, /i is simply the modulation amplitude b averaged over all 
solid angles and the Klein-Nishina cross section. For primary photon energies 
of 10 keV, 100 keV, and 1 MeV the averaged modulation factor equals 0.5, 
0.48, and 0.25, respectively. The angle 0o is the fitted polarization direction 
(in the coordinates of the polarimeter). 

If ot and 0t are the true polarization fraction and polarization direction, 
the probability distribution of the measured values and (po is given by 



3 3,3: 
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In the absence of background, the statistical error on the number of counts 
per bin is given by Poisson statistics: a = y/E. Integration of this distribution 
over 00 gives the probability to measure a polarization degree oq independent 
of dn: 



P{ao) = — ^2 — exp 



4^2 



with Jq the modified Bessel function of order zero. This probability distri- 
bution can be used to determine the minimum fractional polarization which 
can be measured at a certain confidence level. The Minimum Detectable 
Polarization (MDP) is defined as the fractional polarization ai% which will 
be measured in the absence of any true polarization ct = with a chance 
probability of 1%: 

/•oo 

P{a > ai%) = / P(ao) dao = 1% (9) 
Integration gives for a = \ffi: 
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MDP fti — (10) 



As an alternative to the standard method, we evaluate below a Maximum 
Likelihood analysis. The method uses an event list as input: 

L = {(£;,,^,,0,),i = l...no} (11) 

with Ei^ 9i, and (pi being the energy, polar scattering angle, and azimuthal 
scattering angle of the i*^ event, assumed here to be measured with high 
precision. The polarization degree ^^^d polarization direction 0o ^-re de- 
termined by maximizing the Likelihood function: 

no 

A (ao, 0o) = n ^ '^0' (12) 

In practice, it is more convenient to maximize the logarithm of the likelihood 
function In A than to maximize A. The probability of measuring the azimuthal 
scattering angle (pi given the ao, (po, Ei, and 9i is given by the expression: 

P {(Pi; ao, (Po, Ei, Oi) = - (1 + ao^i{Ei, 9i) cos [2(0^ - 0o)]) (13) 

TT 
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with 6i) being the modulation factor for the incident photon energy Ei 
and the polar scattering angle 6i. For the ideal polarimeter, fi is given by Equ. 
(|5]). The maximization of In A gives the maximum likelihood estimates ao,ML 
and </)o,ML of the polarization fraction and polarization direction, respectively. 

We recommend to calculate the statistical significance of the detection 
of a polarized signal by simulating the background and the signal assuming 
the null hypothesis (an unpolarized signal) and calculating the fraction of 
simulated event lists with best-fit polarization degrees exceeding the observed 
one. 

One can estimate the uncertainty on the fit-parameters b y m apping the 



Bayesian posterior probability density function 29| (see also 30|): 



p{ao,(j)o) = (14) 



The function p can be used to plot a region TZa{(3) (or several regions) which 
contains the true value ot and (p^ with a certain probability a: 

TIM = {(«o, 0o) I P {ao, 0o) > (3,aoE [0, 1] , 0o e [0, tt]} (15) 
with /3 is chosen such that 

rfa'p (i0Q p(ao, 0q) = a (16) 



4. Achieved Performance 

We evaluated the performance of the standard analysis and the Maximum 
Likelihood analysis by generating simulated data sets for ax = and the Uq- 
values given in Table [1] For each parameter configuration, 10,000 data sets 
were simulated. We limited the study to primary photon energies of 100 
keV. Individual events were generated by first drawing a polar scattering 
angle from the distribution given by Equ. ([3]). Subsequently, an azimuthal 
scattering angle was generated according to the distribution given by Equ. 
@. The event lists were analyzed with the two analysis methods described 
above. 

Figure [2] shows the distributions of the best-fit polarization fraction for 
the case of no polarization ot = for data sets with no = 10,000 events. The 
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no 


MDP Equ. (dUD 

[%] 


MDP Stand. 
Method [%] 


MDP Max. Lik. 
Method [%] 


1000 


28.3 


28.3 


23.8 


3000 


16.3 


16.2 


13.7 


10000 


8.9 


9.0 


7.4 


50000 


4.0 


3.9 


3.3 



Table 1: Summary of simulation of data sets for a primary photon energy of 100 keV. For 
each number of Compton events, 10,000 data sets were simulated assuming no inherent 
polarization ax ~ 0. For each data set the polarization fraction ao was reconstructed with 
various methods. The ao-distributions were used to determine the Minimum Detectable 
Polarization (MDP) being the value that is larger than 99% of the reconstructed MDP 
values. 

Maximum Likehhood method gives on average smaller best-fit polarization 
values. The MDP is calculated by determining the ao-value which is larger 
than 99% of all fitted values. For no = 10,000 the standard method and the 
Maximum Likelihood method give MDPs of 0.089 and 0.074, respectively. 
Thus, the Maximum Likelihood method results in a reduction of the MDP 
by 21%. It should be noted that the observed MDP for the standard method 
agrees well with the estimate from Equ. ( ITOl) . 

Figure [3] shows the MDPs achieved with the two methods for ng-values 
ranging from 1000 to 50,000 for the Maximum Likelihood method, the stan- 
dard method, and a simpler method. In the case of the "simple method" the 
polarization degree is derived from the maximum number nmax and minimum 
number n-aim of entries in the azimuthal scattering angle histogram according 
to 

ao = . (17) 

^max ~r ^min 

Whereas the standard method uses the contents of all bins in the histogram 
to estimate Oo, the simple method uses only the contents of two bins (the 
bins with the minimum and maximum number of counts). It is thus expected 
that standard method gives more precise estimates of cq. In the absence of 
a true polarization signal, the standard method should thus give on average 
smaller ao- values than the simple method. As a result, the standard method 
leads to a smaller ao-value that can be detected at a certain confidence level 
than the simple method. Figure [3] indeed shows that the simple method 
performs markedly poorer than the standard method. Whereas the MDP 
achieved with the standard method is independent of the number of bins 
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Figure 2: Distribution of the reconstructed polarization fraction ao for a true polarization 
fraction of ax =0. The solid lines show the results for the Maximum Likelihood analysis 
which uses the polar and azimuthal scattering angles of the Compton events; the dashed 
lines show the results for the standard analysis which uses only the azimuthal scattering 
angles (no = 10,000). 



used in the azimuthal scattering histogram (as long as the number of bins is 
not too small to smear out the modulation, and not too large to reduce the 
number of entries per bin to values below 10), the simple method performs 
poorer the more bins are used in this histogram. 

Compared to the simple method and the standard method, the Maximum 
Likelihood leads to a reduction of the MDP by 255% and 21%, respectively, 
independent of uq. As mentioned before, the improvement relative to the 
simple method depends on the number of bins used in the azimuthal scat- 
tering angle histogram. Figure H] shows the distributions of the fitted ao and 
00-values for ot =0.25, 0t = 7r/2, and no = 10,000. The standard analysis 
and the Maximum Likelihood analysis give 68% confidence intervals on oq 
of 29% and 24%, respectively. For the polarization angle 0o the 68% confi- 
dence intervals are 0.06 rad (3.4°) and 0.048 rad (2.7°) for the two respective 
methods. 

5. Summary and Discussion 

Photons of a linearly polarized X-ray beam Compton scatter preferentially 
into the direction perpendicular to the preferred orientation of the electric 
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Figure 3: Minimum Detectable Polarization (MDP) as function of the number of Comp- 
ton events for (from below to above) the Maximum Likelihood method (solid line), the 
standard method (dashed line), and for a simpler method (dotted line) in which the po- 
larization degree is derived from the maximum and minimum number of entries in the 
azimuthal scattering angle histogram. The Maximum Likelihood method leads to an im- 
provement by 21% compared to the standard method. 



field vector. The modulation of the azimuthal scattering angle distribution 
depends on the polar scattering angle 6. The modulation is most pronounced 
for photons with polar scattering angles close to 7r/2 (slightly less for higher 
energies), and weakest for forward-scattered and back-scattered photons. In 
this paper, we have described a Maximum Likelihood analysis which uses 
the polar and azimuthal scattering angles of all detected Compton events 
to estimate the polarization degree and the polarization direction. For an 
ideal Compton polarimeter, the Maximum Likelihood analysis results in a 
reduction of the MDP by 21%. The accuracies ((T68%-values) with which 
non-zero polarization degrees and polarization directions can be measured 
are reduced by similar fractional amounts. Note that the improvement by 
21% results from the properties of Compton scatterings. The improvement is 
somewhat modest, but is consistent with expectations based on Eqs. (jS]) and 
(fTOj) : for events with a certain modulation factor /z the MDP scales inversely 
proportional to /i (Equ. [10]). For the "best events" the modulation factor is 
close to 1; for energies <C mgC^ the average modulation factor is 0.5. The 
decrease of the MDP should thus have a value of between 1 and l//i = 2. 
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Figure 4: Distribution of the reconstructed polarization fraction (left panel) and the 
reconstructed polarization degree (/)o (right panel) for a true polarization fraction of ax = 
0.25 and a true polarization direction (/)t = 7r/2. The solid lines show the results for the 
Maximum Likelihood analysis which uses the polar and azimuthal scattering angles of the 
Compton events; the dashed lines show the results for the standard analysis which uses 
only the azimuthal scattering angles (both panels are for no = 10,000). 



For a real-world Compton polarimeter, a Maximum Likelihood analysis 
could be performed with the Likelihood function given by Eqs. f ll2p and 
f ll3p . For experiments with non- uniform acceptance (e.g. for the ASTRO- 
H soft gamma-ray telescope of the ASTRO-H mission), the PDF of Equ. 
f fT3|) should be replaced with a PDF obtained from Monte Carlo simulations. 
The PDF could then depend on a variety of measured parameters, e.g. the 
estimated energy of the primary event, the location of the first interaction, 
and the polar and azimuthal scattering angles. 

A Maximum Likelihood approach could also be used for the analysis of 
the data from photo-effect polarimeters as the gas polarimeters described in 
0, S, . The approach can be implemented by determining the 
modulation factor as a function of a set Oi of suitable observables for each 
event, and by using the modulation factor /u(Ci) in Equ. (IT^ . One obvious 
choice is to use the reconstructed energy of the events: Oi = {Ei}. Sofiitta 
et al. |32| and Depaola & Longo j33| discuss the pronounced energy depen- 
dence of the modulation factor. Some of the energy dependence stems from 
the atomic physics of the photon-absorption, i.e. the modulation factor for s 
electrons is larger than for p and d electrons and Auger electrons are emitted 
isotropically. Some of the energy dependence stems from the detection tech- 
nique as the properties of the recorded photoelectron tracks depend on the 
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photon energy. One may use additional properties of an event which have 
an impact on the effective modulation factor, e.g., the track length or the 
number of triggered pixels, the ellipticity of the track image when analyzed 
with a second moment analysis, or the estimated error on the reconstructed 
azimuthal scattering angle. The usage of in Equ. f|T3|) will improve the 

sensitivity of the polarimeter as the individual events are weighted according 
to their information content. 
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